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Abstract 

This paper describes, from a deterministic view point, 
an on-line process monitoring method for small signal 
behavior of time invariant stable systems at a given op- 
erating point. The monitoring is accomplished by us- 
ing an m-scqucncc modulated low-pass wavelet as the 
probing signal, and a cyclic reference signal as the cor- 
relator. This method allows to achieve high resolution, 
high signal to noise ratio monitoring, and the acquired 
data can be used for on-line diagnosis. 


1 Introduction 

Health monitoring and diagnosis of industrial processes 
provide situational awareness and feedback informa- 
tion that help optimize operational effectiveness, pre- 
vent catastrophic failures, and decrease maintenance 
costs. Existing approaches to process monitoring in- 
clude model-based, knowledge- based, and multivari- 
able statistical approaches [8]. This paper applies the 
concept of pulse-compression probing to process mon- 
itoring. In this case, a compressible finite duration 
probing signal stimulates a process The response is 
processed by a correlator with a cyclic reference signal 
which compresses the response of the probing signal 
into a narrow source wavelet, called the compressed 
version of the probing signal. Functionally the com- 
pressed version of the probing signal equivalently ex- 
cites the process to yield the probing output from which 
the impulse response that reflects the small signal char- 
acteristics of the process can be extracted. Deviation of 
the impulse response from its normality is an indication 
of abnormality of the process. The acquired informa- 
tion can then be used for further process diagnosis. 

The pulse-compression method has been used in seis- 
mic and acoustic probing applications [2. 3]. This pa- 
per considers the application of the method to on-line 
process monitoring. The principle, the properties of 
the method, the selection of parameters of the prob- 
ing signal, and the implementation issues will be dis- 
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cussed. In addition, its relation to cross-corrclation[5] 
technique in least squares linear system identification 
will be cclincatcd. 


2 Pulse-compression probing method 

Let u(t) = ui(£) + u 2 (t) be the input signal at one 
of the input ports of a stable process, where Ui(t) is 
a prescribed set point or a general reference profile, 
u 2 (£) is a finite duration probing signal, yi(t) be the 
response at one of the output ports of the process to 
u(t) in i he absence of u 2 (t), y 2 (t) be the process re- 
sponse corresponding to the perturbation by u 2 (t), h(t) 
be the impulse response that reflects the small signal 
characteristics of response y to u at the current op- 
erating condition, s(t) be a cyclic reference signal, and 
z(t) = Z{(t) + z 2 {t) be the probing output, where z\ and 
z 2 correspond to y\(t) and y 2 (t), respectively. Let V 
denote the convolution operation, and ‘®’ denote the 
correlation operation. Suppose the probing signal u 2 (t) 
is sufficiently small. Then from Fig.l the following can 
be seen to hold. 

z 2 {t) ■= y 2 {t)®s(t) = h{t)*u 2 {t)®s{t) = h(t)*e{t), (1) 

where e(t) is called an equivalent excitation and is given 
by 

e(t) = u 2 (t) ® s(£). (2) 

If e(£) := £(£), then z 2 (t) = fe(t), and the small sig- 
nal characteristics is extracted. In the following de- 
velopment, we attempt to synthesize e(t) by designing 
signals u 2 (t) and s(t). Note that is fixed by the 

process operation requirement. The signal design prob- 
lem amounts to designing u 2 (t) and s(t) such that the 
following holds. 

y{t)~yi{t), (3) 

or 

rwt.s { } 

rms{yi} 

is minimized, and 

z(t) W z 2 {t) = h(t). (4) 

or 

rms{z 2 } 
rms{zi } 



is maximized. Condition (3) implies that u 2 should be 
constraint in both power and magnitude. The power 
constraint is dictated by the acceptable signal to noise 
ratio in response y(t ), and the magnitude constraint 
is dictated by the requirement to exited the process 
only in its linear region around the current operating 
point. Condition (4) is dictated by the acceptable sig- 
nal to noise ratio in response z(t). Note that the roles 
of signal and noise arc reversed in (3) and in (4). The 
finite duration requirement on the probing signal u(t) 
is consistent with the requirement of a small time delay 
in the probing output that benefits the process moni- 
toring. It is a hard requirement in some applications 
where the duration is limited, for example, by the in- 
terval between the transmitting time and the receiving 
time of the probing signal. 

U ‘2 and s(t) satisfying (3) and (4) can be selected to be 
of the following forms, based on [2. 3]. 


(i) Reference signal s(t) of the correlator. 

The reference signal is the cyclic repetition of an im- 
pulsive m-sequcncc of order n with bit duration to an d 
intensity 1. An example of an impulsive m-sequencc 
m, 6 (t) is shown in Fig. 2 (b) where n = 5, resulting a 
maximum sequence duration of (2° - 1) x to = 31 to. 

(ii) Probing signal u 2 (t). 

u 2 (t) = p{t)*m 6 (t) (5) 

where p(t) is a narrow low pass wavelet of magnitude 
A. Fig. 2 (a) shows two examples of unity magnitude 
wavelets. One is a rectangular pulse and the other 
is a raised cosine pulse. The use of a rectangular 
wavelet will be assumed from this point on for simplic- 
ity. Smoother probing signals such as the m-scquencc 
modulated raised cosine is in general more desirable 
but more complex in signal generation. 


In its original form, a discrete m-scqucncc m(k) of order 
n and magnitude a is a pseudo-random binary sequence 
(PRBS) of maximum length 2 n — 1 with the following 
propcrtics[4, 7]. 


1 


2 n - 1 


2 n — 1 


52 

fc=i 



(6) 


2 n — l 

2 „ ~ 1 52 m (k)m( k + l) 

Ja 2 , 2 k = 0, ±(2" - I), ±2(2 n - 1), - ■ • 
\ — 2^T’ elsewhere 

Note that the second property holds when at least one 
of the two m(A;)’s involved is cyclically rcpcatcd[3]. 

The probing output is given by 

z(t) = y(t) ® s(t) = y x {t) ® s{t) 4- y 2 (t) ® a(t). (8) 


The first term in (8) depends on the fixed input and the 
process. Therefore yi{t) can be assumed to be weakly 
correlated with s(t). Applying (6), the following in- 
equality can be considered valid 


\yi{t) ®s(t ) | < 


max t \yi(t)\p(t - (2 n - l)t 0 ) 

(2 n - l)t 0 


( 9 ) 


Let T = (2 n - l)t Q . The second term of (8) can be 
expressed in a more detailed manner as 


h(t) * m$(t) ® s(t) 

= h(t) * ^{2 n p(t) IT) -p(t) * -jto)}' 

t j 

When n is sufficiently large, we have 

y 2 (t)®s(t)*h(t)* ^- (10) 

to 

If in addition, wavelet p(t) is sufficiently narrow with 
respect to the quickest variations in h(t) the above term 
approximates Ah(t). For the above term to be suffi- 
ciently dominant in z(t ), it must be sufficiently larger 
than max t \yi(t)\/[{2 n -l)to] in (9) which again requires 
sufficien ly large n. 


3 Signal parameters and implementation issues 

Two conditions must be met to arrive at z(t) ss Ah(t). 
These a:*c a sufficiently large n and a sufficiently small 
t Q . In this ease z(t) is called an estimate of h(t), de- 
noted by h(t) . This section will discuss in more rigorous 
terms the meaning of large n and small to- 

Recall that t 0 is the bit duration of the impulsive m- 
scqucnce in use. It is also the width of the equivalent 
excitation e(t) and that of the low-pass wavelet p(t). 
Therefore t 0 determines the level of details of h(t) ob- 
servable from the probing output z(t), or the resolu- 
tion of h(t). In order to observe the full details in h(t), 
to should be determined by using Shannon’s sampling 
theorem [1], which states that if the Fourier transform 
of h(t) is zero outside of (-u^o^’o)^ th cn MO can be 
fully reconstructed based only on its sampled values at 
a sampling rate higher than = 2uq. In practice, h(t) 
is not strictly band limited, and one may settle with a 
sampling rate about 5 times the estimated bandwidth 
ujbw /(27t) of h(t). This rate translates to the condition 

to < 1 /&bw* ( 11 ) 

The probing output repeats itself every T seconds. It 
is obvious that T must be longer than the “memory” 
of h(t) in order to avoid the time domain aliasing. This 
would require that knowledge of how long it takes for 



the slowest significant modes of h(t) to decay to approx- 
imately zero. A reasonable assessment is T = 5t s i ow 
where r siow is the time constant of the slowest mode of 
h(t). Given T and t 0 , n can be determined through 

n > log 2 ( " — b f) (12) 

to 

Therefore, for a specified T, a smaller to (a higher res- 
olution) results in a lager n. In practice, however, 
the above requirement on n is often much less strin- 
gent than that imposed by the signal to noise ratio 
rms{z 2 ] / rms{z\] . For every unit increment in n, this 
signal to noise ratio is doubled. Another direct benefit 
for using a larger n is the reduction in bias as shown in 
(7). Note however, that both a smaller to and a larger 
n lead to an increased computational burden in car- 
rying out the cross-correlation. In addition, from the 
time a probing signal is applied to the time the probing 
output appears in its entirety (one complete h(t), 2 T 
must have gone by. Therefore, whenever the signal to 
noise ratio at z allows, T thus n should be minimized, 
and whenever the resolution on h allows, to should be 
maximized. 

The pulse-compression method for process monitor- 
ing can be easily implemented digitally, though all the 
above discussion was carried out in t.hc continuous do- 
main. The digital data storage solves the problem of 
long delay times normally encountered in using the 
cross-correlation technique on-line for processes with 
large time constants[5]. In this case, the impulses in 
the impulsive m-scqucncc arc replaced by pulses, the 
lowpass wavelet can be replaced by its samples at a 
sampling rate as low as one sample* per bit duration 
(to). The m-scqucncc can be generated by an n-stage 
shift register [4]. The digital signal acquisition and ap- 
plication from and to the process can be accomplished 
by using A/D and D/A converters. 


4 Relation to the cross-correlation technique 

The same pseudo- random excitation is used in both 
the pulse-compression method for process monitoring 
introduced in this paper and in the cross-correlation 
method for impulse response identification for lin- 
ear systems castcd in least squares formulation [6, 5]. 
Therefore the two methods arc in fact equivalent, and 
enjoy the same advantages. Since the methods arc 
nonparametric, little a priory knowledge is required of 
the process model. Processing and memory require- 
ments arc moderate and can be controlled through 
trading off with the monitoring/identification rcsolu- 
tion/accuracy. As long as the process is stable, di- 
vergence is not possible. Though not considered in 
this paper, it has been shown[5. 6| that the cross- 
correlation method has a superb noise rejection ca- 
pability. The pulse-compression method differs from 


the cross-correlation method in the following aspects. 
The pulse compression for process monitoring problem 
is formulated entirely in the continuous domain; the 
probing signal and the reference signal in our pulse- 
compression method are not the same; our method is 
applicable to nonlinear systems; the analysis is carried 
out from a deterministic view; the flexibility of using 
a general lowpass wavelet, (as opposed to rectangular 
only) is introduced; the unncccssity of the periodic rep- 
etition cf the probing signal[7] for less noisy process is 
clarified and therefore reduce the time delay in cross- 
correlating. Use of the monitoring data for diagnosis 
or identification purpose is possiblc[6]. 

For a complete development of the cross-correlation 
method using the m-scqucncc for linear system identi- 
fication, the reader is referred to Hill and McMurtry[5]. 


5 An example 

In this section, the pulse-compression method for 
process monitoring is applied to a simulated electro- 
hydraulic system that is subject to leakage. The system 
exhibits a first order small signal characteristics with 
a bandwidth of approximately 5 rad/scc. In order to 
be able to recognize a small percentage change in leak- 
age coefficient, bit duration t 0 is chosen to be 0.04scc. 
Using the setup of Fig.l, a continuous probing signal 
(u 2 (t)) one tenth of the magnitude of a windowed sinu- 
soidal signal (ui(t)) is applied to the system. u 2 (t) is 
generated through a 15th order impulsive m-scqucncc 
modulated rectangular wavelet. Response y(t) at the 
system output is shown in Fig. 3. The magnitude ratio 
of yi to i/ 2 is about 7. Response z(t) is shown in Fig. 4, 
the magnitude ratio of z 2 to z\ is approximately 40, 
yielding a very accurate estimate of impulse response. 
This is achieved at the expense of a long delay of ap- 
proximately 42(= 2T) minutes before the entire im- 
pulse response is observed. In practice that the noise 
level would not necessitate the pursuit of this accuracy, 
and a shorter delay in acquiring monitoring data is de- 
sirable. 


6 Conclusions 

A pulse-compression method for on-line monitoring of 
the small signal behavior is described. The parame- 
ters in the probing signal can be chosen to balance 
between the speed of processing and the accuracy of 
the acquired monitoring data, which is the impulse re- 
sponse of the process at the current operating condi- 
tion. The paper elucidates the concept of stimulating 
all significant modes of the small signal characteristics 
of a process through the compression of an elongated 
small magnitude poring signal of finite duration into 



a near impulsive equivalent excitation signal. Besides 
the application in process monitoring as discussed in 
this paper, the pulse-compression method has found 
applications in seismic acquisition [2, 3], and in linear 
system identification^, 6, 7] under a different formal- 
ization. Other possible applications include earth crust 
movement monitoring, near surface engineering appli- 
cations, palcontologic study, mechanical malfunction in 
machineries, and biomedical imaging applications. 
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Fig. 2 Lewpass wavelet examples and a probing signal 



Fig.3 Hydraulic system response to sinusoidal 
with/without probing signal 



Fig. 4 Probing output to probing input with/without 
sinusoidal input 


Fig.l A process monitoring schematic 







